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^ ; Abstract 

D ' Rare astrophysical events - cosmological gamma-ray bursts with energies over GeV - are con- 

sidered as an origin of information about some SUSY parameters. The model of generation of the 
. powerful gamma-ray bursts is proposed. According to this model the gamma-ray burst represents 

as a result of the collapse and the total annihilation of the neutralino clump. About 80 % of 
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the clump mass radiates during ~ 100 second at the final stage of annihilation. The annihilation 
spectrum and its characteristic energies are calculated in the framework of Split Higgsino model. 
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I. INTRODUCTION 



Since the direct experimental information about structure of SUSY model and values 
of SUSY particles' masses are absent now we try to extract any possible data from the 
cosmology. From the recent WMAP data the mass density in the Universe is equal 
to its critical value pcrit — 0.54 ■ 10~^GeV/cm^; the structured Dark Matter (DM) and 
the isotropic and homogeneous Dark Energy (DE) contribute to this density ~ 23% and 
~ 73%, correspondingly. As it is known, baryons give ^ 4% of the total Universe's mass. 
In our Galaxy ^ 50% of total mass is the mass of DM which forms the Halo. The Lightest 
Supersymmetric Particles (LSP) - neutralino - are supposed to be the main contributors to 
the cold (and hot) DM. The physics of LSP is discussed intensively in various aspects 01 . 

To observe the neutralino appearance in space or land experiments one can use of the 
different ways. For example, the high energy antiprotons, positrons or neutrino flux due to 
neutralino annihilation in the Dark Halo or in the center of the Sun or the Earth. In our 
work we have concentrated on the specific gamma-radiation J§\ which can be diffuse (due 
to annihilation into fermions and bosons) or monochromatic (due to annihilation into 77 or 
Z'j |3]). To measure characteristics of these specific gamma- fluxes the special experimental 
programs have been elaborated with space telescopes (EGRET, GLAST) 0, @. 

Now we have from the collider experiments only low restriction for the neutralino mass 
> 32.5 GeV. ATLAS and CMS collaborations (LHC) will say something else only after 
2007. But the absolute stability of neutralino (or metastability) cannot be checked in collider 
experiments. Moreover, some details of the standard cosmological scenario, which are based 
on the neutralino as the CDM, should be changed if neutralino's and other superparticles' 
characteristics do not agree with the scenario. Thus we need in combined analysis of the 
LHC and astrophysical data to verify our ideas on cosmological and cosmogonic evolution 



Neutralino in our Galaxy or in other galaxies up to z = 8 10 at the gravity separation 
stage can form massive dense clumps of different sizes 0], 0], j^. Due to the inner gravi- 
tational instability clumps can collapse and intensively annihilate. Our basic supposition is 
that the nature and origin of the rare high-energy GRB can be understood as a result of the 
neutralino clump total annihilation. We have shown that at the final stage of the process the 
powerful gamma-ray burst (GRB) with photon energies over GeV will occur. In the stan- 
dard GRB the energies of photons are only fractions of MeV (see ^3). The specific energy 
and time distribution of this hard 7-radiation can be related with the parameters of super- 
particles (neutralino, in particular) and give an important information on the SUSY model 
details and its compatibility with cosmology. As an interesting example we have considered 
the Split Higgsino scenario fll^ in the framework of famous Split Supersymmetry model 3 



and obtained the values of annihilation time of the neutralino clump, characteristic photon 
energies of the GRB and the annihilation spectrum. 



II. THE COLLAPSE OF THE HOMOGENEOUS NEUTRALINO CLUMP 

Consideration of neutralino as a carrier of DM mass in the Universe generates questions 
about their evolution and about astrophysical phenomenons which prove both the existence 
of neutralino and the reality of the evolution process with neutralino participation. 

According to high resolution cold dark matter simulations in the framework of the hi- 
erarchical principle of structure formation the large virilized halos are formed through the 
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constant merging of smaller halos produced at earlier times. Neutralino at the gravita- 
tional isolation stage 2; = 8 -j- 10 formed dense massive clumps of different scales [q], jl3l |. 
It likes naturally that DM clumps are in the state of dynamical equilibrium as a result of 
compensation of own gravitational forces by forces of tidal interactions with each other. In 
consequence of inner instability the clump can turn into irreversible collapse state. The 
growth of density provides the growth of annihilation rate. So neutralino clouds become 
more and more intense sources of relativistic particles. At final stage of this process there 
is a gamma-ray burst. Specific energy and time distribution of the gamma radiation are 
connected with neutralino parameters and give an important information about structure 
of the MSSM and about its correspondence with cosmology. Therefore our main statement 
consist in the nature of rare powerful cosmological gamma ray bursts is a result of the total 
annihilation of neutralino clumps. 

We have analyzed the qualitative and quantitative descriptions of this process in the 
spherically symmetric collapse model with two assumptions. First, annihilation products 
leave a clump during a time substantially smaller than the time of its macroscopic evolution. 
And second, an annihilating clump is spatially homogeneous and isotropic since more dense 
regions annihilate faster than the rest and heterogeneities are vanishing. 

The model, defined by equations of mass and momentum balance and the equation of 
non-relativistic gravitation theory, is: 



d f /"/ ^ ,r, f ^ i . f ,t . f2{av)annpl ,,(2.1) 

— J p^VidV = - J [p^ViVk + PikjdSk - J p^ViCpdV - J VdisPx^dV - j — VidV , 



A(j) = AnGp^ , 

where Pik is the pressure tensor for neutralino gas. It is absent in consequence of the 
second assumption. Last terms in the balance equations describe the momentary going 
away of annihilation products from the collapsing clump. We have taken into account that 
annihilation and scattering are cross-channels. In the local form we have the simple system: 

dM ?,{av)ann f duV 2GM 



dt 2ttM^ i?3' \dt J R y ^ RoJ' ^^'^^ 

which integrates elementary. Here R is the clump radius; M = Anp^R'^/S is the clump mass; 
Rq = const is the initial radius; < 1 is the constant which fixed by initial velocity of the 
surface clump compression. The results are in the parametric form 
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where M, 



we see M{t) 0, R{t) ^ for t ^ oo. 
The annihilation rate is 



4 ■ lO^^g, p^(*) ~ 10 ^g/cm^ are the initial clump mass and density. As 
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Generally behavior of this function is that a slow growth at small times changes into quick 
evolution around narrow peak after that there is an exponential drop to zero. This behavior 
explains by changing of the parameter hierarchy during the collapse. We use the convenient 
time scale 
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Here to is the time of the maximum. Since we have 
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2ttG 



< 1 



then the annihilation rate doesn't depend on initial conditions and expresses through the 
fundamental constants only. The functions Eann/M{^^) and M/M{^^) of dimensionless time 
1] = {t — to)/ To are shown at the Fig. [T] 
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FIG. 1: Annihilation rate Ea/Mi^^,-^ (right) and relative mass M/Mi^^-^ (left) of neutralino clump as 
a function of dimensionless time r/. 

About 80 % of clump mass annihilates during the time 

tann = ^Tq = . (2.6) 

Effectively the annihilation takes place between ^ = Atq (before maximum) and t^^"* = 2ro 
(after maximum). Such event represents as a powerful gamma-ray burst with total energy 
about clump mass. Kinetic cross section can be parameterized in the form: {(Tv)ann = 
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o? I'^T^Ml^ Mc - some parameter with the dimension of mass. So the decay time of collapsing 
neutralino clump expresses through SUSY parameters Mc, M-^ and gravitational constant 
but it is absolutely independent from initial density and initial velocity of clump compressing. 

Deeper analysis is possible only in the framework of concrete SUSY model. So there exits 
a possibility to extract of the SUSY parameters from observable characteristics of the bursts 



and to find constraints on the SUSY model. In the framework of Split Higgsino scenario [11 



the neutralino mass ~ 3 TeV is estimated from the statement that the Dark Matter in 
the Universe has formed at the high symmetry phase of the cosmological plasma. Neutralino 
of such type in the halo annihilate into W- and Z-bosons and fermions (see Fig. |21). Total 
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FIG. 2: Diagrams of neutralino annihilation in the framework of Split Higgsino scenario, 
kinetic cross section is 

, . _ gt (21 - 40 cos^ Ow + 34 cos^ Ow) , . 

256 7rM2cos4 0VK ' ^ ' 

The clump annihilation time in this model is t^nn — 30 s. 

III. CALCULATION OF CHARACTERISTIC PHOTON ENERGIES AND ANNI- 
HILATION SPECTRUM 

The energy spectrum dN^{E^)/dE^ is determined by the distribution in multiplicity of 
the secondary hadrons resulted from the neutralino annihilation. We use of the fact that 
the total secondary hadrons multiplicity is nearly twice larger than the charged hadrons 
multiplicity [14]. The average multiphcity of secondary charged hadrons {hch){y/s) was 
studied in e~^e~, pp, pp, e^p processes It was established that this quantity is some 

universal function of energy 

fichiVs) = A + Bhiy/s + Chi' ^/s, hch = {rich) (v^/ go) - '^o 

with experimentally fixed parameters 

A = 3.11 ±0.08, = -0.49 ±0.09, C = 0.98 ± 0.02. 

To choice a specific channel it is necessary to fix parameters go 5 ''^o- For neutralino annihi- 
lation go(xx) = 1' '^o(xx) = 0- 

Further it is supposed that the part of charged hadrons n = {rich) / {n-h) doesn't depend 
on energy and has the value k ~ 0.49, which is extracted from Z-peak data jl4| . 
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For characteristic photon energies generation the following processes are most important 
7r° 27, 1]^ — > 37, 77" — > 37r° — >■ 67. In the fermion annihilation channel (Fig. 12)) the total 
hadron multiplicity is described by following logarithmic function with good accuracy: 

(nl^) = K-\A + B ln2M^ + C In^ 2M^) ~ 149 for ~ 3 TeV. (3.1) 

An average energy of neutral pions in neutralino annihilation secondaries is E^^o ~ E^o ~ 
2M;^/ (nf/). Then in the 7r° 27 decay maximal characteristic energy of photon is equal: 

E^(^o^27) - = 20 GeV. (3.2) 

Analogously, for reactions rj^ 37 and 77° 3tt^ 67 energies of photons are: 

E^(^o_3^) ~ ^^0 /3 = 13.5 GeV , ^^^^0^67) - E,^o /6 = 6.7 GeV. (3.3) 

In the boson annihilation channel (Fig. ^ the total hadron multiplicity of W- and Z- 
bosons decays is: 

(n)f ^) ~ 42.9. (3.4) 

Here neutral pion average energy is E^^o ~ Ey^o ~ M^/ (nY^) and maximal characteristic 
photon energies are: 

E^(n"-*2'y) — ~ 35 GeV, 

(3.5) 

E^(^^o_3^) ~ Er,o /3 = 23.3 GeV , ^^(^0^67) ^ E^o /6 = 12 GeV. 

In a wide energy region the multiplicity distribution can be described with a good ac- 
curacy by the Ne gat ive Binomial Distribution (NBD) that depends on energy very weakly 
(logarithmically) |l4l |: 



k~^(^/s) = a + 61n ^/s, 

where n = rich, = "^ch- For various channels the coefficients in the function k~^{^/s) are 
different: 

a^+^_ = -0.064 ± 0.003, fe^+e- = 0.023 ± 0.002; (3.7) 



app/pp = -0.104 ± 0.004, bpp/pp = 0.058 ± 0.001. (3.8) 

Using the NBD it is possible to find the approximate energy distribution of photons. For 
example, in the hadronic channel with Br{h) ~ 0.58 it is supposed that the branching ratios 
for various hadrons Br{i/h) = (jii) / (jih) are nearly constants and equal to the corresponding 
branching ratios extracted from e"^e~-annihilation Each annihilating neutralino pair 
gives a number of neutral pions Br{7T^ / h) ■2nP{n; n, k) with the probability Br{h). These pi- 
ons generate the following distribution of photons An^ ^ 2 ■ Br{h)Br{TT^ / h) ■ 2nP(n; h, k)An 
with energy E^ = M^/2n in the multiplicity interval An, which is connected with energy 
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interval An = M^^E^/2E'^. For 77*^-mesons and boson annihilation channel considerations 
are analogous. Then the number of photons with energy per one annihilation act is: 

^^^7 2M^ ^ [Bri-K^/h) + Br{ri^/h) ■ Br{r]^ ^ 27)] x 



dE. 
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Here Br{WZ) ~ 0.2 is the total branching ratio for neutralino annihilation into W- and 
Z-bosons; charge hadron multiplicities (n^l) and (n^^) are defined in 1)3.11) and ()3.4|1 : param- 
eters kj} = k-\2M^) = 0.4 and fc"^^ = k-\M^) = 0.12 were determined with coefficients 
from ()3.8|) and ()3.7|) . respectively. The spectrum is shown in Fig. El 



dN/dE, 
[l/GeV] 




FIG. 3: Gamma ray annihilation spectrum. 



The particular duration of the GRB and the specific shape of its energy spectrum are the 
characteristic features of considering astrophysical event - the collapse and total annihila- 
tion of the neutralino clump. Calculated values of characteristic photon energies are close 
to the threshold of sensitivity for the satellite detector GLAST. Possibility of neutralino 
annihilation spectrum registration at this apparatus will be clarified when the GLAST will 
be put into operation. 
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IV. CONCLUSION 



We have proposed the new production mechanism of the cosmological gamma-ray bursts. 
In the framework of this mechanism the gamma-ray burst represents as a result of collapse 
and total annihilation of the neutralino clump. The GRB signal when registered should 
contain an information on SUSY parameters and the details of the annihilation process. We 
have calculated the characteristic photon energies and annihilation spectrum for = 3 TeV 
in the framework of Split Higgsino scenario. The particular duration of the GRBs and specific 
shape of its energy spectrum are the characteristic features which can help to separate 
considering astrophysical events from the radiation background. From astrophysical data 
on GRBs with high-energy photons it is possible to extract an information about SUSY 
parameters in different models. It can help to investigate compatibility of the SUSY model 
with cosmology, cosmogony and high-energy collider data. 
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